arXiv: 1508.05919vl [physics.optics] 24 Aug 2015 


Slot-Mode Optomechanical Crystals: A Versatile Platform for Multimode 

Optomechanics 

Karen E. Grutter,^’[^ Marcelo I. Davanco,^ and Kartik Srinivasan^’j^ 

^ Center for Nanoseale Seienee and Teehnology, National Institute 
of Standards and Teehnology, Gaithersburg, MD 20899-6203 USA 
(Dated: August 25, 2015) 

We demonstrate slot-mode optomechanical crystals, a class of device in which photonic and 
phononic crystal nanobeam resonators separated by a narrow slot are coupled through optomechan¬ 
ical interactions. In these geometries, nanobeam pairs are patterned so that a mechanical breathing 
mode is confined at the center of one beam, and a high quality factor (Qo ^ 10^) optical mode is 
confined in the slot between the beams. Here, we produce slot-mode devices in a stoichiometric 
Si 3 N 4 platform, with optical modes in the 980 nm band, coupled to breathing mechanical modes at 
3.4 GHz, 1.8 GHz, and 400 MHz. We exploit the high Si 3 N 4 tensile stress to achieve slot widths 
down to 24 nm, which leads to enhanced optomechanical coupling, sufficient for the observation of 
optomechanical self-oscillations at all studied frequencies. We utilize the slot mode concept to de¬ 
velop multimode optomechanical systems with triple-beam geometries, in which two optical modes 
are coupled to a single mechanical mode, and two mechanical modes are coupled to a single optical 
mode. This concept allows great flexibility in the design of multimode chip-scale optomechanical 
systems with large optomechanical coupling at a wide range of mechanical frequencies. 


I. INTRODUCTION 

Sideband-resolved cavity optomechanical systems have 
recently demonstrated their potential in a wide variety 
of applications, including motion sensing mia, ground 
state cooling EE], and optomechanically-induced trans¬ 
parency For these applications, high efficiency re¬ 

quires large optomechanical coupling strength in addition 
to sideband resolution (mechanical frequency ^ optical 
linewidth). Additional phenomena have been observed in 
multimode cavity optomechanical systems, in which mul¬ 
tiple optical and/or mechanical modes interact, includ¬ 
ing wavelength conversion EU, Raman-ratio thermom¬ 
etry m, energy transfer between mechanical modes m, 
and optomechanical mode mixing m- Phonon pair gen¬ 
eration m, mechanical mode entanglement [12 E], and 
unresolved sideband cooling m have also been theoret¬ 
ically proposed. In all these systems, improved perfor¬ 
mance and broader applicability could be achieved if the 
optical and mechanical modes could be independently 
tailored to a given application. 

The slot-mode optomechanical crystal structure, in 
which the optical and mechanical modes are confined in 
separate but interacting beams (Fig. [^, is one method 
of achieving this flexibility while maintaining large op¬ 
tomechanical coupling strength. Simulations m have 
shown that, in systems in which the optomechanical in¬ 
teraction is dominated by moving boundaries, this geom¬ 
etry can significantly increase the optomechanical cou¬ 
pling strength relative to single nanobeam optomechani¬ 
cal crystals. It also provides the desired design flexibility 
to enable multimode applications such as optomechanical 
wavelength conversion. 
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In this work, we experimentally demonstrate slot-mode 
optomechanical crystals implemented in stoichiometric 
813 X 4 , a material whose broad optical transparency and 
large intrinsic tensile stress make it attractive for many 
applications. In Sec. Ill, we show how this intrinsic stress 
can be exploited to achieve slots with aspect ratios of 
10 : 1 , and in Sec. [IV| we demonstrate how tuning this 
slot width improves device performance in 3.4 GHz band 
devices. Sec. lYl shows how the mechanical mode fre¬ 
quency can be changed while minimally affecting the op¬ 
tical mode, with demonstrations of 1.8 GHz and 400 MHz 
band devices. Finally, in Sec. we extend the slot¬ 
mode optomechanical crystal concept to multimode op¬ 
tomechanical devices, in which two mechanical modes are 
coupled to a single optical mode, and two optical modes 
are coupled to a single mechanical mode. 


II. BASIC DEVICE DESIGN 

A slot-mode optomechanical crystal, shown in Fig. Ef, 
consists of two parallel beams separated by a narrow slot. 
The “optical beam” is a photonic crystal cavity designed 
to confine the optical mode in the slot. The “mechani¬ 
cal beam” is a phononic crystal resonator optimized to 
confine the mechanical breathing mode (Fig. while 
maintaining low optical loss. Both the optical and me¬ 
chanical modes are confined along the z-axis by periodic 
patterning of holes. In the outer mirror region, the lat¬ 
tice spacing is constant, but it varies quadratically in 
the cavity region (Fig. ED- Details on the design of 
this device are outlined in Ref. [16]. The devices in this 
work were designed for optical modes around 980 nm 
and mechanical breathing modes around 3.4 GHz. The 
optical beam is patterned with identical elliptical holes 
(188 nmx330 nm) along its length, while the mechanical 
beam holes have a constant height (370 nm) and widths 
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FIG. 1. (a) Variation of the optomechanical crystal lattice constant along the length of the beams. The period is fixed in the 
mirror regions at the beam ends and varies quadratically in the center cavity region, (b) The slot mode optomechanical crystal 
is formed by parallel optical and mechanical beams that are separated by a narrow slot. The zoomed-in image of the center 
shows the finite element method (FEM) simulated electric field amplitude of the optical slot mode around 980 nm. (c) FEM 
simulation of the breathing mode of the mechanical beam (around 3.4 GHz), (d), (e), and (f): The width of the slot is varied 
in an EEM simulation of the (d) resonant wavelength, (e) optical quality factor (Qo), and (f) optomechanical coupling ^o/(27r). 


that are varied such that the “ribs” between the holes 
align with the elliptical holes in the optical beam. 

There have been several demonstrations of sideband- 
resolved single-nanobeam optomechanical crystals aig 
laHTllIH], in which a GHz frequency mechanical breath¬ 
ing mode is coupled to an optical mode localized by 
the same physical structure. These geometries are dis¬ 
tinguished by the breathing mode’s high frequency (en¬ 
abling sideband resolution), isolation from mechanical 
supports due to the phononic mirrors, and strong inter¬ 
action with the optical mode. Our goal in this work is to 
retain these advantageous features while increasing the 
system’s versatility through the slot mode geometry. 

Optical slot modes have been utilized before to achieve 
large optomechanical coupling in microrings/disks [191 
[20] . bilayer photonic crystal slabs m, and photonic 
crystal zipper cavities [22]. These applications were 
lower frequency (< 150 MHz) than the 3.4 GHz band 
breathing modes in this work, and, thus operated in the 
unresolved-sideband regime (mechanical frequency < op¬ 
tical linewidth). In addition, previous demonstrations 
have not taken full advantage of the flexibility of the slot 
mode architecture, as the mechanical and optical modes 
were supported by the same structural components. Sep¬ 
arating the optical and mechanical modes into two beams 
enables independent design of these modes. This opens 
a wide range of frequency combinations that would be 
difficult to access with a single optomechanical structure. 
The slot-mode structure also opens the possibility for ad¬ 
ditional interactions with both modes, which can be sep¬ 
arately accessed from the beam sides opposite the slot. 
For example, electrodes could be added to the outside of 


the mechanical beam with minimal perturbation of the 
optical mode. Adding more optical or mechanical beams, 
thereby forming more slots, can also increase the device 
functionality by realizing multimode optomechanical sys¬ 
tems, as discussed in Secs. jVlj and m 

III. STRESS TUNING AND DEVICE 
FABRICATION 

In addition to the design of the mechanical and optical 
beams, device parameters are also strongly dependent on 
the width of the slot between the two beams, simulated 
in Fig. [^-f. Given a device with fixed design of the op¬ 
tical and mechanical beams, reducing the slot width red- 
shifts the optical resonance, and reduces Qo somewhat 
(still above 10^). The optomechanical coupling rate go 
increases significantly as the slot width decreases, so the 
slot between the optical and mechanical beams should 
be made as small as possible. Lithographically defining 
small spaces and etching high-aspect-ratio trenches are 
both challenging in fabrication. This can be mitigated 
by taking advantage of the intrinsic film stress of stoi¬ 
chiometric Si 3 N 4 (« 1 GPa). An asymmetric anchoring 
condition in a doubly-clamped beam induces asymmetry 
in the stress, thereby causing it to move laterally. Long, 
thin tethers asymmetrically attached to the ends of par¬ 
allel nanobeams have been used to shrink gaps to as small 
as 40 nm after release [23]. To achieve the same effect, we 
investigated small slits at a beam’s ends (Fig.|^). Finite 
element method (FEM) simulations show that varying 
the width and depth of these slits controls the lateral 
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FIG. 2. (a) FEM simulation of a tensile-stressed beam with stress-tuning slits at the ends, (b) Displacement at beam center 
with respect to slit depth. FEM results (line) are for a beam with the same dimensions as the optical beam of the slot-mode 
device. Error bars on the measured data are due to the uncertainty in the SEM measurements and are one standard deviation 
values, (c) Scanning electron microscope (SEM) images of a released device. Insets show the slot width at the beam end is 
about 70 nm, shrinking to 24 nm at the beam center. 


displacement of the center of the beam (Fig. [^). In the 
slot-mode device, a large initially defined and etched slot 
would be reduced post-release to the desired width by in¬ 
cluding these stress-tuning slits at the ends of the optical 
beam. 

Slot mode optomechanical crystal nanobeams were 
fabricated in 250 nm thick stoichiometric Si 3 N 4 deposited 
via low-pressure chemical-vapor deposition on a bare Si 
substrate (tensile stress 1 GPa). Devices were pat¬ 
terned via electron-beam (E-beam) lithography in a pos¬ 
itive E-beam resist and developed in hexyl acetate at 
7 °C. The pattern was transferred to the Si 3 N 4 using a 
CF 4 /CHF 3 reactive ion etch. Devices were released in a 
45 % KOH solution at 75 °C followed by a dip in a 1:4 
HC 1 :H 20 solution. Finally, the devices were dried on a 
hotplate. 

A scanning-electron microscope (SEM) image of a re¬ 
leased device is shown in Fig. H- The lithographically- 
defined slots were between 80 nm and 120 nm, and, 
with the SEM, we measured stress-tuned slots as small 
as 24 nm at the center, an aspect ratio of about 10:1 
that would be difficult to achieve with lithography alone. 
Fig. graphs the SEM-measured displacements of the 
beam centers with respect to the stress-tuning slit depths. 
The measured trend matches well with the displacements 
predicted in the FEM simulations. 


IV. DEMONSTRATION OF SLOT-MODE 
CONCEPT 

The experimental setup used to characterized the 
Si 3 N 4 slot-mode optomechanical crystals is shown in 
Fig. [3}i, and was previously described in m- All mea¬ 
surements were taken at room temperature and pressure. 
Devices were characterized with a 980 nm external cavity 
tunable diode laser, which was coupled evanescently to 
the devices via a dimpled optical fiber taper waveguide 
(FTW) with a minimum diameter of ~ 1 pm. 

Among the measured devices, a device with a 50 nm 
stress-tuned slot had the highest intrinsic optical quality 


factor Qo at (1.65 ± 0.09) x 10^ (linewidth of 2.0 GHz ± 
0.1 GHz) [24], as shown in Fig. If- QoS up to ^ 2.4 times 
higher have been demonstrated in Si 3 N 4 single-nanobeam 
optomechanical crystals [25], but it is expected that the 
slot mode would have lower QoS because the geometry 
has more scattering sites near the optical mode. In these 
devices, narrower slots generally resulted in lower QoS, 
with 20 nm slot devices having the lowest QoS around 
2.5 X 10^. With further optimization, improving Qo in 
smaller slot designs is feasible. 

We also used optical characterization to more precisely 
determine the effect of the stress tuning. Iterations of 
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FIG. 3. (a) Optical modes are detected by swept-wavelength 
spectroscopy, while mechanical modes are measured when the 
laser is on the blue-detuned shoulder of the optical mode. For 
go calibration, the laser is phase-modulated, (b) Optical res¬ 
onant wavelength of three devices with different stress-tuned 
slot widths, (inset) Optical spectrum and fit of highest mea¬ 
sured Qo among these devices, having designed gap of 50 nm 
and Qo = (1.65 ±0.09) x 10^ [24] (c) Example mechanical 
spectrum, including phase modulator calibration peak. This 
power spectral density plot is referenced to a power of 1 mW 
= 0 dB. Lorentzian fit of thermal noise spectrum is in red. 
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devices were made with the same optical and mechan¬ 
ical design but stress-tuning slits of varying depth, so 
that the only difference among these devices would be 
the final, stress-tuned slot width. An example is shown 
in Fig. Three devices with the same optical and me¬ 
chanical design show a red shift of the optical resonance 
as the designed stress-tuned slot width decreases (the 
stress-tuning slit depth increases). This trend is expected 
from simulation (Fig. W) , and was consistent in 24 of 27 
unique device designs, indicating that varying this stress¬ 
tuning slit depth is a reliable technique for tuning the slot 
width. 



FIG. 4. (a) Mechanical spectra at different input optical 

powers (Pin) for a device with a designed stress-tuned slot 
width of 70 nm, intrinsic Qo = (3.7 ±0.1) x 10^, intrinsic 
= 2380±90 [26], and Qm/(27r) 3.31 GHz. (b) Mechan¬ 

ical spectra at different Pin for a device with 20 nm designed 
stress-tuned slot width, intrinsic Qo — (3.2 ±0.1) x 10^, in¬ 
trinsic Qm = 2400 ± 300, and Qrn/(27r) 3.49 GHz. (c) 

Measured 7m,eff/(27r) of the devices from (a) (blue) and (b) 
(red). Error bars represent the uncertainty in the fit of the me¬ 
chanical spectra to a Lorentzian. Dashed lines show weighted 
linear fits of the subthreshold 7 m,eff/( 27 r). The power spectral 
density plots in (a) and (b) are referenced to a power of 1 mW 
= 0 dB. 

For mechanical mode spectroscopy, the signal was de¬ 
tected with a high-bandwidth (8 GHz) photoreceiver, 
the output of which was sent to a real-time electronic 
spectrum analyzer. Optomechanical characterization re¬ 
quired longer-term stability of the coupling, so the FTW 
was positioned a few hundred nanometers to the side of 
the device and affixed via van der Waals forces to nearby 
protruding parts of the Si 3 N 4 film. The coupling dis¬ 
tance was chosen for a transmission minimum around 
70 %. The blue detuning of the laser further increased 
the measurement stability by enabling access to the ther¬ 
mally self-stable regime [27j so that the laser did not have 
to be externally locked to the cavity. 

We used a calibration signal from a phase modulator to 
measure go in a few devices [28] [29], as shown in Fig. |^, 
where the phase modulator calibration tone is shown 
along with the thermal noise spectrum of the 3.49 GHz 
mechanical breathing mode (quality factor Qm ~ 3900). 
For a device with a designed, stress-tuned slot width 
of 60 nm, we measured ^o/27r = 184 kHz ± 2 kHz, 


where the uncertainty comes from the uncertainty in the 
thermal noise spectrum fit and the measurement of the 
phase modulator = 2.78 V ± 0.01 V (Sec. [^. This 
value matches well with the FEM-simulated go values 
(Fig. [^). Another device, having a designed, stress- 
tuned slot width of 20 nm, had a phase-modulator- 
calibrated ^o/27r = 317 kHz ± 3 kHz, which also aligns 
with FEM simulations and confirms the significant im¬ 
provement in coupling achieved by producing narrow 
slots. We note that such slot mode geometries are of 
particular importance for materials such as Si 3 N 4 , the 
low refractive index (vs. Si or GaAs) of which limits the 
achievable coupling strength in single nanobeam geome¬ 
tries. 

In addition, with the laser blue-detuned (A > 0), op¬ 
tomechanical back-action coherently amplifies the me¬ 
chanical mode, increasing the detected amplitude while 
decreasing the effective mechanical linewidth 7 m,eff- As¬ 
suming only optomechanical damping changes with input 
power, the effective linewidth is related to the optical 
power at the coupling point to the device Pin as follows, 
where k is the intrinsic optical loss rate, /^ex is the exter¬ 
nal coupling rate, Uq is the optical resonant frequency, 
and ftrn is the intrinsic mechanical frequency [301EI]: 

_ 9o ^exPin | hv/2, 

-Im ^ A 2 + {K/2f ( (aTEPW/P 

/t/2 

(A - + {K/2f 

(la) 

~ Tm ± 9o S(/t), /^ex 5 ^oi A, f^m) P’\n 

Thus, the effective mechanical linewidth should change 
linearly with respect to optical power, with the intercept 
indicating the intrinsic mechanical linewidth 7^1 and the 
slope proportional to go. For a blue-detuned laser, this 
slope is negative, and when the optomechanical amplifi¬ 
cation cancels out 7 ^, the device reaches the regime of 
regenerative self-oscillation. The Pin at which this occurs 
is the threshold power. 

We use this relationship to determine the intrinsic 
of these devices by looking at the detected mechanical 
spectrum with respect to power. To compensate for the 
cavity’s power-dependent thermo-optic shift, for each in¬ 
put power, we adjust the laser wavelength to the optimal 
detuning value, which corresponds to the point at which 
the mechanical peak is maximized. We then linearly fit 
the subthreshold 7 m,eff with respect to Pin to find 7 ^ 1 - 
This same procedure is used to compare devices with 
similar optical and mechanical parameters; in this case, 
the slope is an indicator of the relative effective go- 

Fig. 1^ and b show measurements of two such devices 
with similar optical and mechanical performance but dif¬ 
ferent stress tuning. One device, which had a designed 
stress-tuned slot width of 70 nm (stress-tuning slit depth 
of 220 nm), had an intrinsic Qo = (3.7 ±0.1) x 10^ and 
an intrinsic Qm — 2380 ± 90 m- The data correspond- 
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ing to this device are shown in Fig. and the blue 
data in Fig. I§. The other device had a designed stress- 
tuned slot width of 20 nm (stress-tuning slit depth of 
295 nm), an intrinsic Qo = (3.2 ± 0.1) x 10^, and an in¬ 
trinsic Qm = 2400 ±300. The data corresponding to this 
narrower-slot device are shown in Fig. and the red 
data in Fig. H- 

Comparing measurements of these two devices, the me¬ 
chanical mode in the 20 nm slot device is more ampli¬ 
fied than in the 70 nm slot device. For Pin ~ 1.2 mW, 
the detected mechanical peak in the 20 nm slot device 
is ~ 65 dB above the noise floor, while the 70 nm slot 
device’s mechanical peak is only ^ 19 dB above the noise 
floor. The measurements of the optomechanical narrow¬ 
ing of the effective mechanical linewidth (Fig. |^), show 
that the slope of the line for the narrower-slot device (red) 
is much steeper than for the wider-slot device (blue). 
Because they have similar optical and mechanical Qs, 
this suggests that the device with the 20 nm slot has a 
higher effective go. It is also noteworthy that the nar¬ 
rower stress-tuned slot enhances the back-action enough 
that it reaches self-oscillation above a threshold power of 
900 pW. 

Among all the devices measured, devices with more 
aggressive stress tuning (narrower slots) generally had 
steeper linewidth-narrowing slopes, implying higher ef- 
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FIG. 5. (a) FEM simulation of the fundamental lateral flex¬ 
ural beam mode, (b) Measured 3 dB linewidth of the fun¬ 
damental flexural beam mode (blue) and the breathing mode 
(red) as a function of Pin. Error bars represent the uncer¬ 
tainty in the fit of the mechanical spectra to a Lorentzian. 
The fundamental mode self-oscillates at Pin ~ 150 pW, while 
the breathing mode self-oscillates at Pin ~ 900 pW. (c) Side¬ 
bands on the breathing mode (red) and the spectrum of har¬ 
monics of the lower-frequency flexural beam modes (blue) line 
up, indicating a mixing between the two. (inset) The full, 
double-sided spectrum around the self-oscillating breathing 
mode. All power spectral density plots in (c) are referenced 
to a power of 1 mW = 0 dB. 


fective as expected from simulation (Fig. [^) and 
confirmed by the aforementioned phase modulator cali¬ 
bration measurements. The most aggressively tuned de¬ 
vices, with slots designed to be 20 nm, had high enough 
effective gos that all but one of them reached the thresh¬ 
old for self-oscillation within the power range of the laser. 
This indicates that narrowing the slot via stress-tuning is 
an effective way to enhance the optomechanical coupling 
in slot-mode optomechanical crystal devices. 

We note that the measured threshold powers are much 
lower, and the mechanical-linewidth-narrowing slopes 
much steeper, than would be expected with the go values 
obtained from the phase modulator calibration method. 
This suggests other factors in the system are contribut¬ 
ing to the effective optomechanical back-action. These 
could include DC optical gradient forces acting to pull 
the beams closer [22] and interaction of the breathing 
mode with the oscillating flexural beam modes. 

In particular, although we designed these devices for 
optimal coupling to the mechanical breathing mode, and 
focused our measurements on characterizing it, there 
are other mechanical modes that couple to the optical 
mode. Defects in the fabricated device give rise to ad¬ 
ditional breathing-type mechanical modes [32], but the 
most well-coupled modes tend to be the lateral flexural 
beam modes. An FEM simulation of the fundamental 
lateral flexural beam mode (11.4 MHz) of the mechani¬ 
cal beam is shown in Fig. [^. Because it is well-coupled 
to the optical mode and has a much lower frequency 
than the mechanical breathing mode, its threshold power 
for self-oscillation is very low; we measure it to be at 
Pin ~ 150 pW. We also note that, upon detection of the 
optical signal modulated by self-oscillating breathing and 
flexural modes we observed mixing tones as sidebands of 
the breathing mode, as shown in Fig. &• 


V. FLEXIBLE MECHANICAL RESONATOR 
DESIGN 

Separating the mechanical and optical modes into two 
beams in the slot-mode architecture adds flexibility in 
designing these modes compared to a single nanobeam. 
By modifying the design of the mechanical beam, a wide 
range of mechanical frequencies can be accessed with¬ 
out significantly affecting the optical mode. To that end, 
we demonstrate lower-frequency designs around 1.8 GHz 
and 400 MHz. Implementing highly-localized breath¬ 
ing modes in various RF bands (here, the IEEE-defined 
UHF and L) broadens the potential application space. 
One way to change the mechanical frequency is simply 
to change the full width of the mechanical beam while 
keeping the same lattice variation. 

For the 1.8 GHz band design, shown in Fig. [^, the 
mechanical beam width was increased from 700 nm to 
1.55 pm. Measurements of a fabricated device (Fig. [^) 
having a designed, stress-tuned slot width of 80 nm 
are shown in Fig. [^-e. The measured intrinsic Qo = 
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FIG. 6. (a) FEM simulation of the 1.8 GHz band mechanical breathing mode of the 1.55 pim wide mechanical beam, (b) SEM 
image of a fabricated 1.8 GHz band device, (c) Detected mechanical spectra at different ETW input optical powers, (d) At a 
ETW input optical power of 4.7 mW, harmonics of the 1.895 GHz mechanical mode are visible, (e) Measured 7m,eff/(27r) of 
the 1.895 GHz mechanical mode. Error bars represent the uncertainty in the ht of the mechanical spectra to a Lorentzian. The 
dashed line shows the weighted linear fit of the subthreshold 7m,eff/(27r). (f) EEM simulation of the 400 MHz band mechanical 
breathing mode of the 4 pm wide mechanical beam, (g) SEM image of a fabricated 400 MHz band device, (h) Mechanical 
spectra measured at different ETW input optical powers, (i) At a ETW input optical power of ~ 2.6 mW, harmonics of 
the 414 MHz mechanical mode are visible, (j) Measured 7m,eff/(27r) of the 414 MHz mechanical mode. Error bars represent 
the uncertainty in the ht of the mechanical spectra to a Lorentzian. The dashed line shows the weighted linear ht of the 
subthreshold 7m,eff/(27r). The power spectral density plots in (c), (d), (h) and (i) are referenced to a power of 1 mW = 0 dB. 


(1.01 ±0.03) X 10^, and the measured intrinsic Qm = 
2130 ± 50, as derived from the weighted linear ht shown 
in Fig. 1^. These values are comparable to the 3.4 GHz 
band devices. The optomechanical coupling of this de¬ 
vice was strong enough that we observed self-oscillation 
for laser powers above 2 mW (Pin is ^ 20 % of the laser 
power at the FTW input). Above threshold, we also ob¬ 
served harmonics on the breathing mode (Fig. [^). These 
arise from nonlinear modulation of the optical held due to 
the Lorentzian optical mode shape, as reported in other 
systems [331435] . 

For the 400 MHz design, the mechanical beam width 
was increased to 4 pm. At this width, the mechanical 
mode is not well-conhned for the same lattice parame¬ 
ters, but the “ribs” still contribute to the optical con- 
hnement. Thus, we kept the ribs to maintain high Qo^ 
but increased the effective mechanical lattice constants 
by “breaking” two-tbirds of the ribs, as shown in Fig. If. 
Measurements of a fabricated device (Fig. |^) having a 
slot width of 80 nm are shown in Fig. If-j- The measured 
intrinsic Qo = (1.02 ± 0.02) x 10^, and the measured in¬ 
trinsic Qm = 800 ± 300, as derived from the weighted 
linear fit shown in Fig.[^. The Qo is comparable to that 
of the 3.4 GHz band devices, indicating that the “broken- 
rib” geometry minimally perturbs the optical mode. The 
Qm^ however, is much lower than in the 3.4 GHz band 
devices. This is likely due to an increase in air damping 
with decreased frequency [36] and an increase in anchor 


loss from the effective two-thirds decrease in the number 
of lattice periods in the mirror region of the mechanical 
beam. As with the other devices in this work, the me¬ 
chanical spectra (Fig. If) include other, less-well-coupled 
peaks that correspond to additional breathing-type me¬ 
chanical modes from defects in the fabricated device [32| 
or harmonics of lower-frequency flexural modes. The op¬ 
tomechanical coupling of the 414 MHz breathing mode 
of this device was strong enough that it reached self¬ 
oscillation for laser powers above ^1.2 mW. As with 
the 1.8 GHz band device, we observed harmonics on the 
breathing mode above threshold (Fig. [^). 


VI. MULTIMODE OPTOMECHANICAL 
DEVICES 

In addition to increasing flexibility in the available me¬ 
chanical frequencies, the slot-mode device architecture 
enables new functionality in that it is straightforward to 
add another separate optical and/or mechanical mode. 
In this work, we demonstrate two cases: a single optical 
mode simultaneously coupled to two different mechan¬ 
ical beams (“M-O-M”) and a single mechanical mode 
coupled to two different optical modes (“O-M-0”). M- 
0-M devices have a variety of possible applications, with 
theoretical proposals including mechanical mode entan¬ 
glement and phonon pair generation [nun] and ground- 
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FIG. 7. (a) FEM simulation of the optical mode of an M-O-M device designed for coupling to 3.4 GHz band (bottom beam) 
and 1.8 GHz band (top beam) mechanical breathing modes. The optical mode is in both slots simultaneously, (b) SEM image 
of a fabricated M-O-M device, (c) Optical spectrum of M-O-M device. Measurement is in gray, and the Lorentzian fit is in red. 
Measured intrinsic Qo = (1.26 ± 0.02) x 10^ (d) Both mechanical modes measured simultaneously, ETW input power 3 mW. 
Data are in gray and Lorentzian hts are in red. At this optical input power, the 1.93 GHz mode has effective Qm = 3175 di 2, 
and the 3.484 GHz mode has effective Qm = 3350 ± 10, where uncertainty comes from 95 % conhdence interval of ht. (insets) 
EEM eigenmode simulations of corresponding mechanical breathing modes. 


state laser cooling of an unresolved-sideband mechanical 
resonator m- Moreover, recent progress has been made 
in studying M-O-M devices in other platforms experi¬ 
mentally, including recent investigations of Bogoliubov 
mechanical modes [13], as well as systems showing syn¬ 
chronization of mechanical resonators via a travelling op¬ 
tical mode miisH]. An O-M-0 slot-mode device provides 
a new platform for optical frequency conversion, as pro¬ 
posed in Ref. [16]. Unlike in previous demonstrations 
of optomechanically-enabled optical frequency conver¬ 
sion [ 31 ], the O-M-0 device enables quasi-independent 
optical mode selection and independent optimization of 
the coupling into each optical mode. 

In the example M-O-M device of this work, we sur¬ 
round an optical beam with a 1.8 GHz band mechanical 
beam (Sec. 0 and a 3.4 GHz band mechanical beam, 
with 80 nm slots between the beams (Fig.[^). The resul¬ 
tant optical mode is concentrated in both slots simultane- 
ously (Fig.[^). We couple to the optical mode by hover¬ 
ing the FTW a few hundred nanometers above the optical 
beam. The measured intrinsic Qo = (1.26 ± 0.02) x 10^ 
(Fig. [7"), and we simultaneously detect modulation from 
both the 1.8 GHz band and 3.4 GHz band mechanical 
modes (Fig. IT)- For the same input optical power, the 
detected 1.8 GHz band mode has a larger amplitude than 
the 3.4 GHz band mode primarily because a lower fre¬ 
quency mode has a larger thermal noise motional ampli¬ 
tude for the same temperature. With the optical quality 
factor in excess of 10^ (linewidth 2.4 GHz), this device 
is in the range of sideband-resolved operation for both 
the 3.4 GHz and 1.8 GHz band modes, suggesting that 
this device is a candidate for high-frequency Bogoliubov 
mechanical mode studies. We have also measured M- 
0-M devices with 1.8 GHz band and 400 MHz band me¬ 
chanical breathing modes coupled to the same 981.85 nm 
optical mode (Sec, [b]). 

The O-M-0 device demonstrated here comprises a me¬ 
chanical beam with a 1.8 GHz band mechanical breath¬ 
ing mode coupled to an optical beam on each side, with 
80 nm slots between the beams (Fig. It). The top 


nanobeam was made slightly wider, resulting in a red- 
shifted optical resonance. The top and bottom opti¬ 
cal modes were characterized separately by repositioning 
the FTW, and the measured optical spectra are shown 
in Fig. The top mode had a measured intrinsic 

Qo = (1.05 ± 0.02) X 10^ at 973.21 nm, and the bottom 
mode had a measured intrinsic Qo = (1.1 ± 0.1) x 10^ at 
947.34 nm. 

The mechanical breathing mode at ~ 1.835 GHz was 
detected when coupled both to the top and to the bot¬ 
tom optical modes (Fig.[^). (There is another, less-well- 
coupled peak at ~ 1.831 GHz that corresponds to either 
an additional breathing-type mechanical mode from de¬ 
fects in the fabricated device [32] or a harmonic of a 
lower-frequency flexural mode.) We also measured the 
effective mechanical linewidth as a function of power for 
both optical modes (Fig. [^). A weighted linear fit of 
these measurements indicates that the intrinsic Qm as 
measured via each optical mode is in good agreement: 
from the top mode, Qm = 1800 ± 100, and from the bot¬ 
tom mode, Qm = 1800 ± 200. The resonant frequency 
also matches, as shown in Fig. [^, implying that these 
two optical modes are in fact coupled to the same me¬ 
chanical mode. The difference in detected mechanical 
peak heights and the difference in the slopes of mechan¬ 
ical linewidth with respect to optical power stem from 
the fact that the bottom optical mode at 947 nm couples 
more strongly to the mechanical mode. 


VII. DISCUSSION 

We have demonstrated slot-mode optomechanical de¬ 
vices in which the mechanical breathing mode of a pat¬ 
terned nanobeam is coupled to an optical mode that 
is laterally confined by a second patterned nanobeam 
and resides within the slot between the two beams. 
Along with large optomechanical coupling rates in ex¬ 
cess of 300 kHz (as measured via phase-modulator cali¬ 
bration) enabled in part by narrow slot widths that can 
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FIG. 8. (a) SEM image of fabricated O-M-0 device, (b) Separately-measured optical spectra of O-M-0 device. Data are in 
gray, and Lorentzian fits are in red. The 947.34 nm mode (“bottom” beam) has intrinsic Qo = (1.1 ±0.1) x 10^, and the 
973.21 nm mode (“top” beam) has intrinsic Qo = (1.05 ±0.02) x 10^. (insets) FEM simulations of the optical slot modes 
associated with bottom and top optical beams, (c) Mechanical spectra measured at different FTW input optical powers. Top 
spectra were acquired while optically coupled to the top beam, and bottom spectra were acquired while optically coupled to 
the bottom beam, (d) 7m,eff/(27r) as measured via the top optical mode (red) and the bottom optical mode (blue) with respect 
to FTW input power. Dashed lines show weighted linear fits of 7m,eff/(27r). Error bars represent the uncertainty in the ht of 
the mechanical spectra to a Lorentzian. 


be achieved by taking advantage of the tensile film stress 
in Si 3 N 4 , this platform allows for flexible design of the op¬ 
tical and mechanical modes, with mechanical beams tai¬ 
lored to support breathing modes ranging from 400 MHz 
to 3.5 GHz. Moreover, this geometry can naturally be 
extended to multimode systems; we have shown triple¬ 
nanobeam devices with two different mechanical modes 
coupled to a single optical mode, as well as a triple¬ 
nanobeam device in which two different optical modes 
are coupled to a single mechanical mode. 

Future work will focus on the use of these multimode 
geometries in applications such as optical wavelength 
conversion and Bogoliubov mechanical mode formation 


for phonon pair generation. Though some of the current 
devices are already weakly in the sideband-resolved limit 
(/^/27r 2 GHz < f7^/27r ^ 3.4 GHz), additional im¬ 

provements in Qo would enable sideband resolution for 
all of the mechanical frequencies studied. Finally, the im¬ 
plementation of on-chip waveguides will likely be neces¬ 
sary to achieve long-term, stable coupling to multimoded 
systems. 
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Supplementary Material 


This section provides supplementary information for this work. It describes the procedure for determining 
the phase modulator T4-, presents data from an additional triple-nanobeam device, and explores how higher- 
order optical slot modes couple to various mechanical modes. 


Appendix A: Measuring Phase Modulator 14 


In Sec. [IVI 
method [H 


we use the phase modulator calibration 
to measure the optomechanical coupling 
qq of slot-mode optomechanical crystals. The K- of the 
electro-optic phase modulator must be accurately known 
in order to do this calibration. 

To measure 14, we send the 980 nm laser signal 
through the phase modulator, modulated at 3.5 GHz by 
an RF signal generator, and into a scanning Fabry-Perot 
interferometer. The detected signal traces out the car¬ 
rier and phase-modulator-induced sidebands in the opti¬ 
cal signal, and we view them on an oscilloscope, as shown 
in Fig. [ST^ . 
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FIG. SI. (a) An example of the output of the scanning Fabry- 
Perot interferometer for a phase-modulated optical input sig¬ 
nal, as read by an oscilloscope. Spectrum shows the carrier 
peak and the first and second sidebands, (b) Carrier (blue) 
and first sideband (red) peak heights with respect to RF sig¬ 
nal voltage applied to the phase modulator. Points are the 
measured values, with error bars indicating the voltage reso¬ 
lution of the oscilloscope. Lines are the fits of the data. 

Ghanging the power applied by the RF signal genera¬ 
tor to the phase modulator changes the magnitude of the 
carrier and sidebands. The RF power Prf is related to 
the signal voltage Vsig = ^/2ZP^^ where the phase mod¬ 
ulator input impedance Z = 50 11. Knowing this, we can 
graph the peak magnitudes with respect to I4ig, as shown 
in Fig [ST]3 . For a phase modulator, the carrier peak 


magnitude should follow the curve H ( Jq (7 rI4ig/I4))^, 
and the first sideband magnitude should follow the curve 
A{Ji (TrRsig/K-))^, where A scales the amplitude of the 
Bessel functions of the first kind Jq and Ji. We fit data 
from the carrier and first sideband to these functions in 
Fig.[ST]3, and both fits result in 14 = 2.78T0.01 V, where 
the uncertainty comes from the fit and is one standard 
deviation. This value corresponds well with the vendor- 
specified value for the phase modulator. 


Appendix B: Additional M-O-M Device 


In addition to the example M-O-M device in Sec.[^ we 
fabricated and characterized an M-O-M device in which 
a 400 MHz band mechanical beam and a 1.8 GHz band 
mechanical beam surround an optical beam, with 80 nm 
slots between the beams, shown in Fig. [S^ . The resul¬ 
tant optical mode is concentrated in both slots simultane¬ 
ously (Fig. S2 r). We couple to the optical mode by hover¬ 
ing the FTW a few hundred nanometers above the optical 
beam. The measured intrinsic Qo was (5.80 ± 0.06) x 10^ 
(Fig. |S2|3), and we simultaneously detect modulation of 
the transmitted optical signal from both the 400 MHz 
band and 1.8 GHz band mechanical modes (Fig. [S^). 
For the same input optical power, the detected 400 MHz 
band mode has a larger amplitude than the 1.8 GHz band 
mode primarily because a lower frequency mode has a 
larger thermal noise motional amplitude for the same 
temperature. As with the other devices in this work, the 
mechanical spectrum (Fig.|S2|i) includes other, less-well- 
coupled peaks that correspond to additional breathing- 
type mechanical modes from defects in the fabricated de¬ 
vice [3] or harmonics of lower-frequency flexural modes. 


Appendix C: Coupling to Higher-Order Slot Modes 

These slot-mode optomechanical crystals confine mul¬ 
tiple optical modes in the slot in addition to the fun¬ 
damental mode for which they were designed. We ob¬ 
served some of these higher-order modes, as shown in 
Fig. [S3^ and b. Because they are distributed more widely 
along the slot, these modes couple less strongly to the 
highly-localized breathing mode and more strongly to 
other mechanical modes in the device, such as higher- 
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FIG. S2. (a) FEM simulation of the optical mode of an M-O-M device designed for coupling to 1.8 GHz band (bottom beam) 
and 400 MHz band (top beam) mechanical breathing modes. The optical mode is in both slots simultaneously, (b) SEM image 
of a fabricated M-O-M device, (c) Optical spectrum of M-O-M device. Measurement is in gray, and the Lorentzian fit is in red. 
Measured intrinsic Qo = (5.80 ± 0.06) x 10^ (d) Both mechanical modes measured simultaneously, ETW input power 2 mW. 
Data are in gray and Lorentzian fits are in red. At this optical input power, 400 MHz mode has effective Qm = 1030 ± 20, and 
1.927 GHz mode has effective Qm = 1450 ± 20, where uncertainty comes from 95 % confidence interval of fit. (insets) EEM 
eigenmode simulations of corresponding mechanical breathing modes. 
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EIG. S3, (a) Eundamental optical slot mode at ~ 949 nm with 
intrinsic Qo = (3.96 ± 0.09) x 10^. (b) Higher-order optical 
slot mode in the same device at ~ 971 nm with intrinsic Qo = 
(3.57 ±0.08) X 10^. (c) Mechanical spectra measured while 
coupled to the fundamental optical mode (red) and higher- 
order optical mode (black). The optical power input to the 
fiber taper waveguide was ^2.2 mW at the 949 nm mode and 
~ 3.5 mW at the 971 nm mode. This power spectral density 
plot is referenced to a power of 1 mW = 0 dB. 


order breathing-type mechanical modes arising from fab¬ 
rication defects [3]. 


In the example of Fig. a device with a designed, 
stress-tuned slot width of 20 nm has an optical mode at 
949 nm with an intrinsic optical Q = (3.96 ± 0.09) x 10^ 
as well as an optical mode at ~ 971 nm with an intrin¬ 
sic optical Q = (3.57 ± 0.08) x 10^. When pumped at 
the 949 nm fundamental optical mode, it self-oscillates 
at the 3.52 GHz mechanical breathing mode, and we 
begin to see sidebands due to mixing with the modu¬ 
lation from the low-frequency flexural beam modes, as 
described in Sec. 4 of the main text. However, even with 
50 % more optical power, the mechanical breathing mode 
does not self-oscillate when pumped at the 971 nm optical 
mode. In addition, the mechanical spectrum reveals an¬ 
other peak at ~ 3.48 GHz with about the same optome¬ 
chanical coupling to the 971 nm mode as the ~ 3.52 GHz 
mode, suggesting that this higher-order optical mode is 
also coupling to some higher-order, less-well-confined me¬ 
chanical mode. 
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